function x = num_mutation(x,pm)

if rand < pm
        x=x + (1-2*rand)* 10e-3 ;
    %x= x + [(1-2*rand) * 10e-3  (1-2*rand) * 10e-3 ];
end
